Quantum loop models and the non-abelian toric code 
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I define quantum loop models whose degrees of freedom are Ising spins on the square lattice as 
in the toric code, but where the excitations should have non-abelian statistics. The inner product 
is topological, allowing a direct implementation of the anyonic fusion matrix on the lattice. It also 
makes deconfined anyons possible for a variety of values of the weight per loop d in the ground state. 
For d = y/2, a gapped non-abelian topological phase can occur with only four-spin interactions. 



Non-abelian anyons have been the subject of intense 
study recently, especially because of their potential ap- 
plication to topological quantum computation It is 
now well understood how abelian fractionalizcd excita- 
tions can occur in relatively simple spin systems, e.g. 
the "toric code" Q, and quantum dimcr models An 
essential property for having deconfined anyons in such 
models is that the ground state can be expressed as a 
superposition of states comprised of closed loops of all 
lengths. Anyons are attached to each other by segments 
of loop, so that their braiding is non-trivial when far 
apart. There thus has been considerable effort to find 
generalizations which have fractionalized excitations with 
non-abelian statistics in "quantum loop models" 0]. 

The Hilbert space of a quantum loop model is spanned 
by loop configurations on some lattice. A Hamiltonian 
of Rokhsar-Kivelson type [j| is chosen so that when the 
ground state is written as a sum over different loop 
configurations, the coefficient for each loop configuration 
is the Boltzmann weight of some classical loop model. 
Here I study the "completely-packed" loop model, where 
every link of the lattice is covered by self- and mutually- 
avoiding loops. There is therefore a quantum two-state 
system at every vertex, corresponding to the two ways 
possible for the loops to avoid each other [3, @] , shown 
in fig. Q] Each loop configuration C with n c loops re- 
ceives weight g?™ c in the ground state. The correspond- 
ing classical loop model is the Q = d 2 -stdXe Potts model 
at its self-dual point: the loops surround clusters in the 
Fortuin-Kasteleyn expansion [7(. 

When defining a quantum loop model, the inner prod- 
uct as well as the Hamiltonian must be specified. The 
simplest inner product is to make each loop configura- 
tion an orthonormal basis element [J. This, however, is 
undesirable because of the u d = \/2" barrier. Correla- 
tors in the ground state of a quantum-mechanical model 
arc weighted by l^l 2 . When each loop configuration is 
orthonormal to the others, 
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i.e. each loop gets a weight d 2 . The classical loop parti- 
tion function ( V E'| V I / ) is dominated by "short loops" when 
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FIG. 1: The two-state quantum system at each vertex 



the weight per loop d 2 is larger than 2 |9| . This short-loop 
phase is not critical, so loops of arbitrarily-long length do 
not appear in the ground state of the quantum model. 
This means that anyons are confined when d > \/2. 

There is another reason why the simplest inner prod- 
uct is undesirable [|| . Non-abelian anyons have loop seg- 
ments attached to them. Consider two such states with 
anyons at the same locations. One then can "glue" the 
dangling loop ends of the two together, so that the com- 
bined configuration consists entirely of closed loops. To 
obtain a topological theory in three-dimensional space- 
time, each loop formed by this gluing must contribute a 
factor of d to the topological inner product, just as loops 
in the ground state do. For example, let \rf) and |x) each 
have four anyons in the same places, but let the loop 
ends be connected in different ways. Computing the in- 
ner product is easily done via the schematic pictures in 
figure El giving (rj\ri) = (x\x) = d 2 , while = d. 





FIG. 2: Gluing to find the topological inner product 

Loop configurations are not orthonormal with the topo- 
logical inner product. 

In this paper I define the completely-packed quantum 
loop model on the square lattice with a local and topo- 
logical inner product. With this inner product, the parti- 
tion function ( x f , | x f ( ) of the corresponding classical model 
is no longer given by fl}, but by a different classical loop 
model. There is good evidence that the loops in the new 
classical model are critical for values of d larger than s/2, 
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so it is possible to kill two birds with one inner product. 

The topological inner product in the completely- 
packed quantum loop model is found by considering the 
entire system to be a single vertex, so that |0) and |0) 
are the only two states in the system. Each has four 
ends, and the topological inner product is computed by 
gluing ends at the same site together. Gluing (0| to |0) 
gives two loops, as in the left of figure [2 Each loop 
contributes a factor d, so the 1/d in front in figure Q] nor- 
malizes (0|0) = (0|0) = 1. Gluing (0| to |0) gives a single 
loop, as in the right of figure The demand of topolog- 
ical invariance therefore requires that at every vertex 



(0|0) = <0|0> = 1, (0|0) = 1/d. 
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This inner product is positive definite when |<2| > 1; note 
that d can be negative. 

This choice of inner product is very natural. As illus- 
trated in figure [31 define 

ID = -^=(d\0)-\0)) (3) 
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This yields (0|1) = (0|1) = and (1|1) = (1|1) = 1. 
There are therefore two natural orthonormal bases for 
the Hilbert space at each site, one with basis elements 
|0),|1), and the other with basis elements |0),|1). The 
unitary transformation relating them is 



\T) ~ F ID 
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This is precisely the desired fusion matrix for non-abclian 
anyons in the quantum loop model, or equivalcntly, the 
fusion matrix of the conformal field theory SU(2)k with 
d = 2cos(7r/(fc + 2)) for k integer. 

+ - J r--f 

FIG. 3: |1) (up to overall normalization) in terms of loops; 
1 1} has each term rotated by 90 degrees. 



The reason why this is exactly what we want goes to 
the very heart of the quantum loop model. Excitations 
in the quantum loop model are obtained by cutting a 
loop, so that two quasiparticles are attached by a zero- 
energy strand. The weighting of d per loop in the ground 
state is necessary for the excitations to have consistent 
non-abelian braiding and fusing (in mathematical lan- 
guage, these are the conditions for a conformal field the- 
ory/modular tensor category) Namely, in the SU(2)k 
theories considered here, these consistency conditions re- 
quire that each type of anyon corresponds to a repre- 
sentation of the quantum-group algebra U q (sl(2)). One 



doesn't need to know much about this algebra other than 
that it has spin-s representations like ordinary su(2) for 
s < k/2, and that their tensor product behaves similarly 
that of sw(2), once one takes this truncation into account. 
In particular for k > 2, two spin-1/2 representations fuse 
to either spin or spin 1. With these rules, like with 
su(2), there are two linearly-independent ways that four 
spin-1/2 representations can fuse to spin 0. 

The anyonic structure is beautifully realized in this 
quantum loop model. Since anyons are attached to 
strands, the strand must fuse as if they have spin 1/2. 
Two anyons connected by a strand fuse to the identity: 
the combination has trivial statistics. Demanding that 
the strands form closed loops in the ground state thus re- 
quires that the four strands at each vertex on the square 
lattice fuse to the identity. Labeling the four strands at 
a vertex (or equivalcntly, the four spin-1/2 representa- 
tions) by a, b, c, d, the state |0) corresponds to requiring 
that the pair a, b fuses to the identity (spin 0). Since all 
four strands fuse to the identity, the pair c, d must then 
also fuse to the identity (0x0 = 0). Likewise, the state 
|0) corresponds to the pair 6, c and the pair a, d each fus- 
ing to the identity. As if these were su(2) spins, these two 
states are not orthogonal. Instead, the state orthogonal 
to |0) corresponds to the pair a, b and the pair c, d each 
fusing to spin 1, which is denoted as the state |1). The ba- 
sis |0), |1) is defined analogously using the pairs a, d and 
fe, c. The change-of-basis matrix F in the fusion algebra 
is the same as that obtained in J5]). Thus the topological 
inner product indeed requires that the strands behave as 
the spin-1/2 representation of SU(2)k- 

This Hilbert space and ground state are related to 
those of "string-net" models [l(|, where non-abelian 
anyons arise by fine-tuning the Hamiltonian to reflect the 
fusion matrix of the desired topological field theory. For 
the SU{2)k string-net model, the states allowed on each 
link of the honeycomb lattice are labeled by the represen- 
tations of spin j = 0, 1/2, . . . , k/2. Analogously to [H and 
here, the ground state is comprised solely of states where 
the representations on the links touching each vertex all 
fuse to the identity. The Hilbert space of this paper is 
obtained by restricting the SU(2)k string-net model so 
that all the links in two of the three directions of the 
honeycomb lattice have j = 1/2, because requiring that 
every trivalent vertex fuse to the identity means that the 
states on the remaining links can take the states j = or 
1. They correspond to the states |0) and |1) on the square 
lattice simply by compressing each such link to a point. 
The |0),|1) basis at a vertex corresponds to stretching 
out each square vertex in the orthogonal direction. 

Now that the Hilbert space and the inner product 
have been specified, I construct a Rokhsar-Kivelson-type 
Hamiltonian (i.e. "rokk" the Hilbert space) whose ground 
state is a sum over loop configurations with weight d per 
loop. This Hamiltonian is a sum over projection oper- 
ators which annihilate the ground state, and all eigen- 
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states with non-zero eigenvalues arc orthogonal to the 
ground state. The off-diagonal "flip" part acts on pla- 
quettes with all four links belonging to the same loop, as 
illustrated in figured On the plane, one can connect any 
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FIG. 4: The flip term F P~P P~ 



configuration to any other by doing flips like this and its 
rotations by 90 degrees 0, [q] . 

The projection operators which implement the flip and 
enforce the weight d per loop in the ground state therefore 
have off-diagonal terms flipping between |0) and |0). In 
the |0), |0) basis, these operators are 
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where the subscript indicates the configuration with 
larger weight in the ground state. Operators in the |0) , |0) 
basis have a hat and those in the |0), |1) basis have a tilde. 
The projectors in the orthonormal bases are found from 
(j6|) by the non-unitary change of basis 
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In the orthonormal |0), |1) basis the flips are therefore 
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A slightly confusing fact is that the operator which 
projects onto |0) in the |0), |0) basis is not diagonal when 
transformed to the |0), |1) basis. Rather, 
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(P oc ^ jjj is normalized so that P Q is a projector.) 

The off-diagonal pieces are a consequence of (0|0) 7^ 0: 
P should be understood as projecting onto a state or- 
thogonal to |0), i.e. in the |0), |1) basis 

P = fp q f = p t . 

Likewise, the operator P~ projects onto a state orthogo- 
nal to |0), i.e. |1): 



v T Kv . = ( J J) =p 1 







The Hamiltonian acts on a two-state "spin" system 
|0), |1) at each site of the lattice, each term involving 
four spins around a plaquette. Then H = ^ p F v , with 

T v = J 7 ° o-^o °o M)° 0*^0 °o°o°o^o 
where the first projector in each term acts on the lower- 
left spin on the plaquette, the second on the upper-left, 
and so on clockwise. By construction, the sum over all 
configurations with weight d per loop is annihilated by H . 
Since this Hamiltonian is the sum of projection operators, 
all eigenvalues must be greater than equal to zero. When 
space is a plane, repeated applications of H connect all 
the configurations, so the ground state is unique. 

On the torus, however, the ground state is not unique, 
because J- p will not change the number of loops which 
are wrapped around a cycle of the torus - it only cre- 
ates and annihilates loops locally. To get a finite number 
of ground states on the torus, one must therefore add 
another term to the Hamiltonian. A local term which 
breaks this degeneracy while still preserving the weight 
d per loop in the ground state is possible only when 
d = 2 cos(7rj'/(fc 4 2)) with j and k + 2 coprime inte- 
gers This term is known as the Jones- Wenzl pro- 
jector, and in the algebraic picture, it projects onto the 
state of spin (k 4 l)/2, whose corresponding anyon is 
not part of the spectrum. For d = \2, fusing two spin- 
1 anyons gives only the identity sector (only the spin- 
1/2 anyon is non-abelian when k = 2). Therefore, fus- 
ing three spin-1 anyons here cannot give the identity: 
1 ® (1 ® 1) = 1 (8) = 1. Since the topological characteris- 
tics of the ground state should reflect this fusion algebra, 
loop configurations involving the three spin-1 strands in 
figure [5] should be forbidden from the ground state. 




FIG. 5: The Jones- Wenzl projector TqTqPqTq for d = ±\/2 



The projector onto the state pictured in figure[5]is nice 
in the square-lattice model here, since it involves only the 
spins around a single plaquette. It involves only projec- 
tors given above, because the flip operators J- and T~ 

introduced in (J6]) project onto |1) and |1) respectively: |1) 
is the eigenstate of the projection operator Tq with non- 
vanishing eigenvalue. Therefore, the Jones- Wenzl projec- 
tor for a plaquette in the |0), |1) basis at k = 2 is 

where the first projector in each term acts on the lower- 
left spin as before. Adding this to the Hamiltonian en- 
sures that the combination of loops illustrated in figure 
flSJ) is orthogonal to the ground state. 
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The full Hamiltonian for d = \[2, with nine ground 
states on the torus, is then given by H = J2 P + Jp\ ■ 
It can be written out entirely in terms of Pauli matrices 
a x and a z acting at each site, with each term involving 
four Pauli matrices acting around a plaquette. Applying 
the F matrix ([5| at d = V2 means that changing be- 
tween |0), |1) and |0), |1) bases amounts to interchanging 
<j x and a z in the Hamiltonian. This resembles the toric 
code Q if one treats each spin as living on the links of 
another square lattice with unit cell of twice the area as 
the original lattice. Then the Hamiltonian divides into 
terms acting on the four links around each site, and the 
four links around plaquette of this new lattice. 

One way of obtaining non-abelian anyons in the spec- 
trum is by allowing empty links (defects) into the model. 
However, expanding the Hilbert space is probably un- 
necessary. Because Tq projects onto spin 1, the net spin 
of a plaquette of spins not annihilated by T v is also 1. 
For k > 2 spin-1 anyons have non-abelian statistics, but 
even for k = 2, these excitations can braid non-trivially 
because of the two spin-1/2 strands attached to each. 
It is not yet proven whether or not such excitations are 
gapped (the proof of gaplessness in Q does not apply be- 
cause of the different inner product) . It seems very likely 
that once the Jones- Wenzl projector is included, they will 
be gapped. Each J p can be added to the Hamiltonian 
with any coefficient without changing the ground state, 
so for an excited state not to be gapped, it would need 
to be annihilated by all J v like the ground state. 

With the topological inner product, it is possible to 
crack the d = \/2 barrier. To have deconfined anyons, the 
corresponding classical loop model must be critical 0, @] • 
The classical partition function is no longer simply (JTJ) , 
but instead is given by a sum over two completely-packed 
loop configurations C and £' on the same lattice: 



d n cd n cl> n x 



(7) 



where n x is the number of vertices at which C and C 
differ. With the topological inner product, A = 1/d, 
while (fTJ) from the simple inner product is recovered when 
A — > 0. Since taking the inner product of a given configu- 
ration with itself always gives an even power of d, chang- 
ing d — > —d is equivalent to instead flipping A — > —A. 
(Note that in the orthonormal bases the Hamiltonian 
depends only on d 2 .) Since the partition function of 
the Q-state Potts model can be expanded in terms of 
completely-packed loops with weight y/Q each, de- 
scribes two eP-state Potts models coupled by a self-dual 
perturbation. When A = ±1, the models are decoupled. 
For A — * one obtains a single Q = d 4 state Potts model 
at its self-dual point, which indeed is critical only for 
d< V2. 

When d = \[2, the partition function ([7]) is the same as 
that for the Ashkin- Teller model, i.e. two coupled Ising 
models. For any A > 0, including A = 1/d = 1/V2, the 



model remains critical [9(. Along this critical line, the 
energy operator (odd under duality) is of dimension x = 
2/g, where g = 8/tt sin -1 [1/2 + 1/(2 + 2y/2X)]. However, 
understanding the fractal properties of the loops as a 
function of A seems to be an open problem. 

As opposed to the classical loop model with partition 
function |T]), the model with can be critical even for 
\d\ = 2cos(7r/(fc + 2)) > s/2 if A is negative. By us- 
ing level-rank duality of the loop representation of the 
BMW algebra [ll[, it is shown in [12] that there is a 
critical point when A = A c , where A c = —y/2sm(n(k — 
2)/[4(fc + 2)]). Moreover, numerical evidence strongly 
suggests that for 1/A C < A < A c , the classical loop model 
has a critical phase [ljj]- When k = 6, A c = — 1/d. Thus 
for k < 6, the classical loop model at A = —1/d falls into 
this phase, deconfining anyons in the quantum model. 

This quantum loop model can be generalized by relax- 
ing the requirement that all links be in the same repre- 
sentation [1, or by studying different lattices. The 
completely-packed model on the square lattice discussed 
here simplifies in several very nice ways, e.g. at d = \[2 
the Jones- Wenzl projector is no more complicated than 
the flip term. On the Kagome lattice at d = y/2 it should 
be even simpler, only involving three-spin interactions. It 
would be very interesting to find other lattices and other 
representations with such elegant properties. 
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